clc
clear all

% Define the x and y values
x = [2 3 4 5 7];
y = [0.61 1 2.12 0.43 0.015];  % Adjusted y values to match the length of x
y = y*1.6;  % Adjusted y values to match the length of x

% Generate a finer set of x values for smoother plotting
x_fine = linspace(min(x), max(x), 100);

% Perform pchip interpolation to get smoother y values
y_smooth = pchip(x, y, x_fine);

% Create the plot
figure;
plot(x_fine, y_smooth, 'k--', 'LineWidth', 5);  % Black dashed line for the smooth curve
hold on;
plot(x, y, 'ko', 'MarkerFaceColor', 'w', 'LineWidth', 2, MarkerSize=30);  % Original data points with white-filled circles

% Set the figure to be square
axis square;



% Set axis limits to better match the example plot
xlim([1.5 7.5]);
ylim([-0.2 3.7]);

% Display grid
% grid on;

% Set specific plot properties to match the uploaded image
set(gca,'fontsize', 35, 'LineWidth', 3);  % Adjust font size and axis line width
set(gca, 'xtick', [2:1:7]);
% set(gca, 'ytick', [2:1:7]);
% Add labels and title
xlabel('Bulk pH', 'Interpreter', 'tex',FontSize = 45);  % Using LaTeX interpreter for the sigma symbol
ylabel('Fe/ SPM Molecule', 'Interpreter', 'tex',FontSize = 45);  % Greek Gamma and units in brackets